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1.  INTRODUCTION 


Carbon-carbon  (C-C)  composites  have  established  themselves  as  impor¬ 
tant  high-temperature  materials  for  spacecraft  applications.  Proposed  new 
aerospace  designs  such  as  those  in  the  National  Aerospace  Plane  (NASP)  and 
in  advanced  turbine  engines  also  utilize  C-C  composites  extensively. 

Some  C-C  composites  are  graphitic  materials.  Such  materials,  includ¬ 
ing  carbon-fiber-reinforced  car bon -matrix  composites,  are  generally  manu¬ 
factured  at  temperatures  exceeding  2400®C.  Glassy-carbon-matrix  C-C 
composites  also  are  frequently  heat-treated  to  relatively  high  tempera¬ 
tures,  above  1700“C,  to  achieve  their  maximum  mechanical  properties.^  When 
such  advanced  composite  materials  are  used  for  applications  requiring  high 
thermal  stability,  thermal  strain  becomes  an  important  engineering 
property . 

The  bulk  thermal  strain  of  C-C  composites  is  ideally  describeo  in 
terms  of  the  properties  of  the  constituent  materials.  Improved  composites 
can  then  be  designed  without  the  laborious  trial -and-error  of  manufacturing 
test  composites.  However,  C-C  composites  are  complex,  and  this  type  of 
predictive  design  is  difficult  to  apply.  Thermal  strain,  the  reversible 
dimensional  change  associated  with  heating,  is  an  approximately  linear 

p 

expansion  for  most  solids  but  not  for  graphite.  Graphite's  basal  plane 
has  a  minimum  strain  near  ^00®C,  and  there  is  a  large  crystallographic 
thermal  expansion  anisotropy  between  the  basal  plane  and  the  c-axis.  These 
factors  make  precise  thermoelastic  modeling  of  a  composite  difficult. 

Also,  data  are  lacking  on  the  pyrolysis  and  high-temperature  behavior  of 
precursors  and  matrix  carbons.  During  the  manufacture  of  a  C-C  composite, 
the  carbon  matrix  precursor  is  pyrolyzed  and  carbonized  to  convert  it  to 
carbon  or  graphite.  In  the  process,  the  matrix  densifies  and  a  complex 
crack  network  and  stress  field  form  in  the  matrix  and  at  the  matrix- 
filament  interfaces,^  making  it  difficult  to  form  an  a  priori  estimate  of 
the  final  thermal  strain  of  a  given  composite  design.  Further  investiga¬ 
tion  of  the  matrix  pyrolysis  behavior  and  its  correlation  with  bulk  thermal 
strain  is  therefore  required. 


5 


A  problem  arises  in  laboratory  investigations  of  such  phenomena  over 
the  2000 -to- 3000 ®C  range:  The  upper  temperature  limit  of  most  commercially 
available  high-temperature  dilatometers  is  only  about  2000*0.  We  therefore 
constructed  and  tested  a  high-temperature  graphite  dilatometer  that  can 
accurately  measure  thermal  strain  phenomena  at  temperatures  as  high  as 
3000*0.  Before  this  apparatus  was  constructed,  high-temperature  dilatomet- 
ric  and  rheological  measurements  could  not  be  carried  out  at  The  Aerospace 
Corporation,  and  such  facilities  are  generally  not  available  at  commercial 
testing  laborator’es. 

This  report  describes  the  dilatometer  and  is  organized  as  follows. 

The  primary  components  of  the  dilatometer  system  are  described  in  Section  II. 
The  theoretical  and  practical  aspects  of  dilatometer  calibration  and  the 
method  of  data  reduction  are  discussed  in  Section  III.  Finally,  initial 
experiments  to  measure  thermal  strain  on  several  carbonaceous  materials 
with  the  dilatometer  are  sunmarized  in  Section  IV. 
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II.  DESCRIPTION  OF  APPARATUS  AND  SAFETY  CONSIDERATIONS 


The  dilatometer,  diagrammed  in  Fig.  1,  consists  of  a  sample  tube/push- 
rod  and  a  cage  assembly  that  are  fabricated  from  a  single  billet  of  poly¬ 
crystalline  graphite  to  ensure  uniformity  of  thermal  expansion.  The  sample 
cage  can  be  removed  from  the  sample  tube  so  that  other  cages  can  be  fabri¬ 
cated  and  installed  to  acconmodate  samples  of  different  sizes.  The  sample 
pushrod  extends  coaxially  through  the  sample  tube  to  make  contact  with  the 
sample.  The  mass  of  the  pushrod  (-10  g)  serves  to  keep  it  in  contact  with 
the  sample  during  the  test.  A  linear  variable  differential  transformer 
(LVDT),  whose  core  is  fixed  to  the  pushrod,  measures  the  displacement  of 
the  pushrod  (as  the  sample  expands  upon  heating)  and  sends  an  electrical 
signal  to  a  data  acquisition  system.  The  heating  elements,  which  are  also 
made  from  graphite,  are  located  coaxially  around  the  sample  cage. 

The  furnace  (manufactured  by  Astro  Division  of  Thermal  Technologies 
Inc.)  consists,  from  the  heating  elements  outward,  of  a  graphite  retort 
around  which  is  packed  thermal  insulation  (lampblack  and  graphite  felt), 
encased  in  an  airtight,  water-cooled  aluminum  housing.  The  interior  of  the 
aluminum  housing  is  purged  with  an  inert  atmosphere  during  the  test  to 
prevent  oxidation  of  the  insulation  and  retort,  heating  elements,  sample 
holder  assembly,  and  sample.  The  furnace  case  and  the  LVDT  are  cooled  by 
flowing  water  from  the  laboratory  cooling-tower  water  supply. 

The  cooling  water  has  an  emergency  city-water  (tap  water)  backup.  If 
the  cooling- tower  water  flow  is  interrupted,  e.g.,  by  a  power  outage,  two 
valves  can  be  manually  opened  to  provide  city-water  cooling  to  the  furnace 
case.  The  cooling-tower  water  lines  are  monitored  by  a  flow  sensor  that 
sounds  an  alarm  (100-dB  intensity)  and  shuts  off  the  furnace  power  when  the 
cooling-water  flow  rate  drops  below  normal.  Following  such  an  interrup¬ 
tion,  the  time  allowed  at  a  given  temperature  before  the  aluminum  furnace 
housing  fails  frcrni  overheating  has  been  estimated  by  heat  flow  calculations 
(see  Fig.  2).  At  a  maximum  operating  temperature  of  2500“C,  5  min  would  be 
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Fig.  1.  Cutaway  view  of  the  dilatometer,  showing  LVDT  and  sample  tube/ 

pushrod  assembly  in  place.  The  top  closure  of  the  furnace,  which 
supports  the  LVDT,  and  the  entire  aluminum  furnace  case  are  water- 
cooled. 


8 


Fig.  2.  (a)  Heat  flow  out  of  furnace  case  as  measured  from  change  in 

the  temperature  of  the  cooling  water  flowing  at  a  known  rate, 

Q  =  me  AT:  m  is  the  cooling-water  flow  rate,  c  is  the  heat 
capacity  of  water,  AT  is  the  temperature  rise  of  the  cooling 
water,  and  Q  is  the  heat  flow  into  the  furnace  shell  ffoc  the 
interior,  (b)  Approximate  times  for  the  exterior  of  the 
furnace  to  reach  melting  temperature.  Note  also  that  these 
calculations  assume  constant  heat-flow  rate,  i.e.,  that  the 
interior  of  the  furnace  stays  at  constant  temperature  after 
water  interruption,  which  it  does  not.  The  actual  time  to 
failure  is,  however,  regarded  as  50?  shorter  than  indicated  in 
order  to  allow  a  safety  margin. 
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available  for  the  operator  to  shut  down  the  power  so  that  no  potential 
safety  problem  is  presented.  The  possibility  of  a  cooling-water  failure 
dictates  that  furnace  operation  be  carried  out  only  during  the  normal 
working  shift,  limiting  rur  length  to  about  8  hr.  A  "tag-team”  system  is 
used  to  ensure  the  continuous  presence  of  a  qualified  individual  in  the 
laboratory  during  that  shift.  Figure  3  is  a  diagram  of  the  cooling  system 
with  the  emergency  bypass  system  noted. 

Sample  dilation  and  temperature  data  are  recorded  simultaneously  by  a 
multitasking  laboratory  computer,  a  Hewlett-Packard  Integral  PC.  The  data 
from  completed  runs  are  sent  as  ASCII  files  for  analysis  to  a  VAX  n/785 
mini-mainframe  computer  via  the  Intelligent  Data  Switch  (IDS).  The  Inter¬ 
active  Data  Language  (IDL)  prograiosing  language,  used  for  data  analysis  on 
the  VAX  computer,  greatly  reduces  programming  and  data  analysis  time  while 
retaining  the  inherent  flexibility  of  user -designed  data  analysis  software. 
The  BASIC  language  code  for  data  acquisition  is  included  in  Appendix  A,  and 
the  IDL  code  for  data  analysis  is  in  Appendix  B. 

The  temperature  of  the  furnace  and  that  of  the  test  sample  are 
measured  by  two  different  systems;  For  temperatures  below  1()00®C,  a 
Chrome 1-Alumel  thermocouple  is  positioned  in  the  sample  cage,  Just  below 
the  sample.  For  the  900-to-3000®C  range,  optical  pyrometry  is  used  and 
requires  two  separate  pyrometers:  a  total -emission-type  pyrometer,  which 
covers  the  range  1000  to  1550®C  and  measures  the  total  flux  of  radiation 
from  a  specific  area  of  the  sample;  and  a  so-called  two-color  pyrometer, 
which  reads  from  1550  to  3000®C  and  measures  the  intensity  ratio  of  two 
adjacent  near- infrared  wavelengths.  For  both  pyrometers,  the  sample's 
emissivity  is  assumed  to  approximate  that  of  a  blackbody.  For  maximum 
accuracy,  a  small  cavity-radiator  serves  as  a  target  for  the  pyrometer 
measurements.  The  cavity-radiator  was  used  to  reduce  the  undesirable 
reflectivity  that  the  system  may  have  on  the  sample  so  that  better  accuracy 
can  be  achieved.  One  side  of  the  radiator  is  located  within  1  mm  of  the 
sample  during  the  run  to  maintain  near -iso thermal  conditions.  As  the 
furnace  temperature  passes  1000®C  during  a  typical  run,  the  thermocouple  is 
withdrawn  manually.  The  temperature-measurement  function  automatically 
passes  from  the  thermocouple  to  the  optical  pyrometers  at  about  950‘‘C. 
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Fig.  3 
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.  Diagram  of  the  cooling-tower  water  system  for  the 

dilatometer  furnace.  If  the  system  fails,  valves  A  and 
B  are  closed  and  valves  D  and  C  are  opened  ( in  that 
order)  to  provide  a  city-water  backup. 
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The  selection  of  the  particular  inert  purge  gas,  helium,  for  the 
dilatometer  is  important  for  very-high-temperature  applications  with  an 
electric  field  present.  At  temperatures  above  about  25^0°C,  argon,  the 
most  conmonly  used  of  the  noble  gases,  begins  to  ionize  and  conduct 
electricity.  The  location  of  both  electrodes  at  the  same  end  of  the 
heating  element  in  this  particular  furnace  design  results  in  a  large 
electric  potential  gradient  there  (-50  V/cm).  When  attempts  were  made  to 
heat  the  furnace  above  25^0‘’C  in  an  Ar  atmosphere,  a  spontaneous  short 
developed  across  the  heating  elements  and  the  furnace  would  cool.  The 
short  carried  about  300  A  at  10  V  and  was  attributed  to  thermally  assisted 
ionization  of  the  Ar.  When  helium  was  substituted  in  the  furnace 
atmosphere,  the  short  disappeared  completely. 

The  useful  upper  temperature  range  of  the  dilatometer  is  dictated 
primarily  by  the  lifetime  of  the  polycrystalline  graphite  heating  element. 
At  3000®C,  the  lifetime  of  the  element  is  less  than  30  min;  below  this 
temperature,  the  life  of  the  element  increases  substantially.  Wachi  and 
Gilmartin^  found  that  at  high  temperatures  (>2900®C),  the  carbon  binder 
phase  in  a  polycrystslline  graphite  begins  to  evaporate  preferentially, 
generating  gaps  between  the  graphite  grains;  the  graphite  grains  are  then 
(apparently)  ejected  from  the  degrading  material.  Electric  arcs  that  may 
form  in  the  resulting  gaps  in  the  dilatometer  heating-element  material 
might  also  accelerate  the  erosion  process.  The  service  life  of  the  heating 
element  when  the  furnace  is  run  to  3000®C  is  about  three  weeks,  and  that  of 
the  furnace  retort  is  about  three  months.  Upper  temperatures  are  therefore 
usually  limited  to  2850®C. 
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III.  CALIBRATION  OF  THE  INSTRUMENT 


A.  THEORETICAL 

The  accuracy  of  the  single-pushrod  diiatometer  design  is  limited  by 
the  thermal  expansion  of  the  materials  out  of  which  the  dilatcxoeter  is 
fabricated.  Vitreous  silica  is  the  preferred  material  for  high-precision 
dilatometers  because,  for  many  applications,  its  coefficient  of  thermal 
expansion  (CTE)  contributes  negligibly  to  the  overall  expansion  of  the 
sample /diiatometer  system.  However,  vitreous  silica  dilatometers  must  be 
corrected  for  high-accuracy  measurements  on  materials  whose  thermal  expan- 
sions  are  close  to  that  of  silica  glass  itself  (about  5  x  10  '  °C  ). 

Silica  cannot  be  used  in  this  diiatometer,  however,  because  the  useful 
maximum  temperature  for  silica  is  only  about  1100®C.  The  diiatometer 
discussed  here,  therefore,  is  constructed  of  a  fine-grained  poly crystalline 
graphite  (SPEER,  Inc.)  to  permit  high  temperatures  to  be  investigated.  The 
average  CTE  of  such  graphites  ranges  from  about  3  to  7  *  10"^  from 
room  temperature  (RT)  to  3000®C,  which  is  about  an  order  of  magnitude 
greater  than  that  of  silica.  Consequently,  the  diiatometer  must  be  care¬ 
fully  calibrated  to  remove  the  apparent  CTTE  contribution  of  the  probes  in 
the  sample  pushrod/tube  assembly. 

The  procedure  we  used  for  correcting  the  diiatometer  data  arises  from 
the  mechanics  of  the  diiatometer  system  and  is  generalized  for  single- 
pushrod  dilatometers;  it  will  now  be  described.  Thermal  strain  occurs  at 
specific  locations  in  the  diiatometer,  as  indicated  in  Fig.  Aa.  For 
simplicity,  the  figure  includes  only  the  primary  structural  parts  of  the 
diiatometer.  The  three-point-bend  fixture,  which  replaces  the  dilation 
sample  in  the  hot  zone  for  flexural  studies,  is  shown  in  Fig.  4b.  The  same 
correction  procedure  can  be  used  in  the  flexure  mode  if  an  inflexible 
standard  sample,  such  as  pyrolytic  graphite,  is  used  for  the  calibration. 

The  sample  and  the  part  of  the  pushrod/tube  assembly  that  are  in  the 
hot  zone  are  heated  during  a  test  and  undergo  thermal  strain.  The  total 
strain  is  measured  by  the  LVDT  at  C.  The  contribution  of  the  "overlapping" 
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(b) 


Fig.  4.  (a)  Schematic  of  single-pushrod  dilatometer. 

(b)  Detail  of  three-point-bend  fixture.  (F  is  force.) 
During  use,  the  fixture  is  positioned  in  sample  zone  A 
of  (a). 


portion  of  the  probes  (in  region  B)  to  the  total  measured  strain  is  zero  if 
the  materials  that  compose  them  are  identical  in  thermal  expansion,  and  if 
the  temperature  gradients  along  the  tube  and  the  pushrod  are  the  same. 
However,  the  thermal  strain  of  the  tube  and  pushrod  do  not  cancel  in  the 
region  where  the  sample  is  located  (D);  the  thermal  strains  from  both  the 
sample  (A)  and  the  sample  tube  in  region  D  contribute  simultaneously  to  the 
net  thermal  strain  measured  at  the  LVDT,  Therefore,  the  LVDT  actually 
measures  the  thermal  strain  difference  between  the  sample  length,  A,  and 
the  corresponding  length  of  the  sample  tube,  D.  The  arrangement  of  the 
probes  requires  that  the  thermal  strain  of  the  length  of  the  sample  tube  in 
region  D  be  added  algebraically  to  the  total  measured  thermal  strain, 
because  the  strain  of  this  part  of  the  tube  acts  to  move  the  center  of  mass 
of  the  sample  downward,  away  from  the  LVDT.  The  thermal  strain  character¬ 
istic  of  the  sample  tube,  at  least  in  the  region  overlapped  by  the  sample, 
must  first  be  precisely  known  before  this  can  be  done. 

The  characteristic  uniform  thermal  strain  of  the  sample  tube  should  be 
ascertainable  by  measuring  the  thermal  strain  of  a  well-characterized 
reference  material.  However,  it  was  found  that  the  thermal  strain  of  the 
graphite  sample  tube  could  not  be  normalized  with  respect  to  change  in  the 
sample  length  (see  Fig.  7).  Therefore,  measurements  of  a  single  reference 
sample  of  fixed  length  were  not  adequate  to  characterize  the  tube  for  all 
sample  lengths.  (The  correction  for  this  effect  will  be  discussed  in  the 
next  section.)  The  difficulty  appears  to  be  due  to  the  large  and  unequal 
thermal  gradient  in  the  sample  tube  and  pushrod  at  the  edge  of  the  hot 
zone,  and  the  fact  that  different  lengths  of  the  pushrod  enter  the  hot  zone 
for  different  sample  lengths. 

To  bypass  this  problem  without  extensively  characterizing  the  baseline 
associated  with  various  sample  lengths,  the  sample  length  was  standardized 
at  1  cm.  The  measured  strain  at  the  LVDT,  then,  can  be  thought  of  as  the 
sample  strain  plus  the  dilatometer 's  instrumental  response: 

NET  , 

“sample  ■  “sample  “machine  ^ 
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which  can  be  rewritten 


“sample  ■  '“sample  ■  “STD^  “STD  “MACHINE  ' 

where  is  the  known  expansion  of  a  standard  reference  material.  A 

sample  of  the  standard  reference  material  can  be  run  in  the  dilatcwneter, 
yielding,  by  Eq.  ( 1) , 


NET 

“STD  '  “STD  “machine 


(3) 


NET 

The  net  strain  of  the  standard  (“gyp)  csui  be  subtracted  from  that  of  the 
unknown,  Eq.  (2),  to  yield 


NET 

“corrected  *  “sample  “  “STD 


(it) 


Since  was  assumed  to  be  known,  it  can  be  added  (as  previously 
mentioned)  to  this  result  to  yield  the  true  thermal  strain  of  the  unknown 
sample. 

The  material  for  a  standard  sample  should  be  selected  because  it  is 
either  precisely  characterized  or  it  possesses  a  low  enough  CTE  that  any 
imprecision  in  the  characterization  of  its  CTE  is  of  little  significance. 
Practical  candidates  are  silica  (fused  quartz  glass),  Zerodur-type  glass 
ceramics  (5  *  10"®  ®C"^),  and  titanium  silicates  (<  5  *  10“*^  ®C"^). 
Partially  crystallized  glasses,  such  as  Zerodur,  are  limited  to  about 
600®C.  Silica  can  be  subjected  to  almost  1100®C  and  can  be  made  more 
refractory  and  lower  in  CTE  by  the  addition  of  up  to  titania  (ULE).® 
The  melting  ranges  of  titanium  silicate  glasses  have  not  been  clearly 
identified  in  our  laboratory,  but  such  glasses  appear  to  withstand 
temperatures  of  about  1H00®C  without  melting,  although  devitrification 
occurs  at  the  sample  surface.  In  general,  the  above  materials  offer  such 
low  expansions  that  even  an  error  as  large  as  50f  in  their  CTE  data  repre¬ 
sents  an  order-of-magnitude-lower  error  than  is  expected  from  the  graphite 
dilatometer.  (Silica-based,  single-pushrod  dilatometers  of  this  type  are 
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practically  limited  to  measurements  of  the  order  of  10"^  m  due  to 
mechanical  errors  in  the  apparatus.^)  National  Bureau  of  Standards  (NBS, 
now  referred  to  as  NIST,  National  Institute  of  Standards  and  Technology) 
standard  reference  material  (SRM)  738  (sapphire),  which  is  no  longer 
available  from  NIST,  can  also  be  used  since  it  is  very  accurately  cali¬ 
brated,  but  its  temperature  range  is  limited  to  about  ISSO^C. 

B.  PRACTICAL  ASPECTS  OF  DILATOMETER  CALIBRATION 

Polycrystalline  graphite  has  a  very  high  use  temperature,  and  conse¬ 
quently  it  has  been  employed  as  a  standard  material,  in  other  work,  after 
careful  optical  calibration.^  However,  pyrolytic  graphite  (PG)  has  a  lower 
CTE  that  is  considered  to  make  it  a  better  standard  material  because  errors 
in  estimating  its  CTE  would  be  correspondingly  low.  The  theoretical 
(crystallographic)  expansion  of  graphite  parallel  to  the  basal  plane  is 
about  that  of  silica  (1  x  10"^  ®C"^),  a  fact  that  favors  its  selection; 
however,  owing  to  the  crystalline  imperfections  inherent  in  PG,  an  expan¬ 
sion  about  two  times  higher  than  this  is  usually  observed. 

For  the  above  reasons,  a  standard  PG  sample  was  used  in  this  dilatom- 
eter  and  was  developed  in-house  as  follows.  A  7-cro-long  sample  of  PG  was 
pre-anneal ed  at  2750‘’C  for  2  hr  in  a  graphite  tube  furnace,  then  reheated; 
its  thermal  strain  was  measured  optically  (above  SSO'C)  using  traveling 
telemicroscopes  that  were  sighted  on  fiducial  marks  in  the  ends  of  the 
rod.  A  1-cm  piece  was  excised  frcxn  this  sample  for  use  as  a  standard.  It 
was  then  calibrated  to  750*’C  in  a  DuPont  943  thermomechanical  analyzer 
(silica  dilatoraeter) ,  which  was  referenced  to  NBS  silica  (SRM  739).  These 
thermal  strain  data  were  then  consolidated  to  yield  the  thermal  strain 
behavior  of  the  sample  over  the  entire  temperature  range  of  the  dilatom- 
eter.  A  polynomial  curve  was  then  fit  to  the  data  for  use  in  subsequent 
data  analysis.  Figure  5  plots  the  data  from  this  standard  and  the  polynom¬ 
ial  fit.  Published  data  for  the  thermal  strain  of  ATJ-S  poly crystal line 

Q 

graphite  are  plotted  in  curve  (b).  Fig.  5,  for  reference. 
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Fig.  5.  Thermal  strain  data  from  (a)  optical  calibration  of  the  in-house- 
developed  pyrolytic  graphite  standard,  (b)  published  data  for 
ATJ-S  polycrystalline  graf^lte,®  and  (c)  theoretical  basal-plane 
expansion  of  graphite.  The  dashed  curve  (d)  is  the  corrected  PC 
data  from  the  in-house  standard. 


18 


Measurements  of  the  NBS  738  sapphire  sample  in  the  dilatometer  after 
this  calibration  procedure  indicated  the  need  for  applying  the  following 
correction  to  the  optical  calibration  data  for  the  PG  standard.  The  CTE 
data  for  a  sample  of  NBS  738  sapphire,  together  with  the  published  data,^ 
are  plotted  in  Fig.  6.  None  of  the  materials  involved,  in  either  the  stan¬ 
dards  or  the  equipment,  undergoes  phase  changes  during  heating,  and  all  are 
considered  homogeneous  and  isotropic.  Therefore,  it  was  assumed  that  these 
materials  display  expansion  characteristics  described  by  a  function  that  is 
continuous,  strictly  increasing,  and  has  first  and  second  derivatives  that 
are  also  continuous  and  strictly  increasing  (i.e.,  the  curve  is  very 
"smooth,"  having  no  inflections).  This  is  a  good  assumption  since  graphite 
is  nonpolymorphic  at  the  tempe'.'atures  and  pressure  used,  and  pre-heat 
treatment  to  high  temperatures  renders  it  stable  at  lower  tempera¬ 
tures.  In  order  to  make  these  smoothness  assumptions,  the  dilatometer 

had  been  preconditioned  by  heating  to  3200®C  for  0,5  hr  prior  to 
calibration. 

The  deviation  of  the  measured  sapphire  standard  data  from  the 
published  data  increased  very  smoothly  (see  deviation  data  in  Fig.  6a), 
thereby  indicating  that  a  linear  error  was  present  in  the  optical  calibra¬ 
tion  data  of  the  PG  standard.  The  indicated  correction  for  this  deviation 
was  made  to  the  optical  calibration  data  by  extrapolating  the  deviation  to 
3000®C  and  subtracting  the  resultant  curve  from  the  PG-standard  data 
(Figs.  5  and  6), 

To  observe  the  dilatometer 's  fine-scale  instrumental  response,  samples 
of  ATJ-S  polycrystalline  graphite  were  measured  in  the  dilatometer. 

Because  the  probes  are  also  made  of  poly crystalline  graphite,  the  net 
thermal  strain  measured  by  the  LVDT  was  predicted  by  Eq.  (1)  to  be  c  zero 
except  for  such  fine-scale  nonlinearities.  Uncorrected  [i.e.,  Eq.  (1)] 
expansion  curves  that  illustrate  this  behavior  were  obtained  from  three 
ATJ-S  graphite  samples  of  different  lengths  from  the  same  block  of  graphite 
and  are  plotted  in  Fig.  7.  Because  the  thermal  expansion  of  the  graphite 
in  the  samples  and  probes  is  assumed  to  be  very  smooth,  the  inflected  shape 
of  the  curve  is  attributed  to  interaction  of  the  thermomechanical  and 
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Fig.  6.  (a)  Initial  dilatooeter  measurement,  following  calibration,  on 

NIST  sapphire  standard  reference  material  SRM  738  (solid  line). 
The  published  data  for  this  material  are  displayed  pointwise  as 
open  circles.  The  deviation  between  measured  and  published  data 
is  also  shown  (dashed  line).  The  indicated  error  is  in  the 
optical  calibration  data  (>  700®C). 

(b)  Two  successive  runs  on  SRM  738  after  subtracting  deviation 
data  (extrapolated  to  3000*’C)  frcmi  graphite -standard  calibration 
curve.  (See  also  Fig.  5.)  Published  data  for  SRM  738  are  also 
shown,  as  open  circles. 
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Fig.  7.  Thermal  expansion  data  for  ATJ-S  poly crystalline  graphite  samples 
of  different  lengths.  Since  the  ATJ-S  graphite  has  almost  the 
same  CTE  as  the  sample  tube/pushrod  assembly,  the  uncorrected  data 
have  a  small  net  slope.  The  inflections  present  in  the  curves  are 
attributed  to  the  instrumental  response  of  the  dilatometer,  since 
the  samples  and  the  sample  probes  undergo  no  phase  changes  during 
heating.  The  dilatometer's  response  to  changes  in  sample  length 
is  nonlinear  and  not  aonotonically  varying  with  respect  to  sample 
length.  The  length  of  sample  for  each  curve  is  indicated. 
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electronic  systems  in  the  dilatometer.  Although  it  varies  with  sample 
length  in  a  complex  manner,  the  fine-scale  instrumental  response  can  be 
considered  constant  if  the  sample  length  varies  only  by  the  amount  typi¬ 
cally  associated  with  the  thermal  strain  of  high-temperature  materials 
(<20  X  10"^  [For  very  large  strains,  such  as  those  associated  with 

pyrolyzing  carbon  precursor  materials  (see  Fig.  12),  fine-scale  nonlinear¬ 
ities  can  be  neglected.]  Such  nonlinearities  are  included  in  the  instru¬ 
mental  response  term  [Eq.  (3)]  by  using  a  higher-order  polynomial  function 
to  fit  the  data  prior  to  subtraction,  and  a  seventh-order  polynomial  fit 
was  found  to  be  adequate.  The  room-temperature  sample  length  also  must  be 
held  constant;  we  therefore  standardized  the  room- temperature  sample  length 
at  1  cm.  Based  on  the  data  shown  in  Figs.  5  and  6,  the  resolution  of  the 
system  is  better  than  ±5  wm,  giving  a  typical  sensitivity  of  about 

0.5  X  10*^  “C"''. 

The  above  rationale  and  data-correction  procedures  are  also  applied  to 
analyze  data  from  dilatometer  runs  with  bar  samples  in  the  flexure  mode, 
using  the  arrangement  shown  in  Fig.  4b. 
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IV.  RESULTS 


In  this  section  we  triefly  describe  the  results  from  early  runs  of  the 
dilatometer  to  illustrate  the  utility  and  flexibility  of  the  system.  The 
discussion  is  also  intended  to  inform  about  the  general  thermal  strain 
behavior  of  carbonaceous  materials. 

Uncorrected  data,  representing  the  net  thermal  strain  of  the  dilatcMD- 
eter-sample  system  [Eq.  (1)],  are  presented  in  Fig.  8  for  non-heat-treated 
PC  and  the  PG  standard  sample.  The  slope  is  negative  because  the  sample 
data  are  uncorrected.  The  progressive  stabilization  of  the  PG  (which  is 
manufactured  at  about  24(X)°C)  and  the  expansion  of  the  PG  standard  (which 
had  been  pre-annealed  at  2750°C)  are  indicated.  The  magnitude  of  the 
sharp,  irreversible  increase  in  the  basal-plane  dimension  of  the  non-heat- 
treated  sample  diminishes  progressively  as  crystalline  defects  are  annealed 
out  and  the  degree  of  graphitic  order  of  the  sample  approaches  a  limiting 
value^®*^^  with  increasing  time  at  high  temperature.  The  GTE  of  the 
pyrolytic  graphite  in  the  basal  plane  changes  little  as  a  result  of  the 
high-temperature  heat  treatment,  indicating  good  stability  of  the  sample's 
CTE  with  respect  to  use  of  the  standard  at  high  temperature. 

Typical  expansion  data  taken  from  various  carbon  materials  are  shown 
in  Figs.  9  and  10.  The  true  (corrected)  data  t^ken  both  perpendicular  and 
parallel  to  the  basal  plane  of  a  sample  of  pyrolytic  graphite,  Fig.  9, 
exhibit  the  large  thermal  strain  anisotropy  characteristic  of  graphite,  as 
well  as  the  irreversible  strains  associated  with  further-graphitizing  heat 
treatment,  past  the  original  temperature  of  manufacture. "nie  data  in 
Fig.  10  are  for  representative  samples  of  the  three  principal  weave  geom¬ 
etries  of  C-C  composites:  one-,  two-,  and  three-dimensional  weaves.  They 
are  very  similar  in  CTE  because  their  thermal  strains  are  dcmiinated  by 
their  constituent  carbon  filaments,  which  vary  in  CTE  by  only  about  a 
factor  of  2  from  that  of  basal -plane  graphite. 
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Fig.  8.  Uncorrected  thermal  strain  data  from  pyrolytic  gr^hite.  Curves 
(a)  and  (b)  are  from  successive  runs  on  as-received  samples; 
curves  (c)  and  (d)  are  from  the  in-house-calibrated  standard.  The 
CTE  of  the  pyrolytic  graphite  (excluding  its  initial  high-tempera¬ 
ture  pyrolysis  expansion)  is  relatively  constant  even  after  the 
large  (irreversible)  strains  that  accompany  the  initial  heatings 
above  the  temperature  of  manufacture  (about  2400”C). 
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Fig,  9.  True  thermal  strain  response  of  pyrolytic  graphite  in  (a)  basal- 
plane-normal  direction  (c-axis)  for  the  second  run,  and  (b)  and 
(c),  parallel  to  the  basal  plane  (a-axis)  for  the  first  and  second 
runs,  respectively.  The  irreversible  shrinkages  and  expansions 
occurring  in  the  c-  and  a-axes  above  2500*C  are  associated  with 
further  graphitization.  Curve  (c)  indicates  that  the  sample 
becomes  stable  with  respect  to  thermal  strain  after  sufficient 
heat  treatment. 
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Fig.  10.  Thermal  strain  of  several  carbon-carbon  composite  configurations: 
(a)  a  Cartesian  223  weave  three-dimensional  composite  manufac¬ 
tured  by  General  Electric;  (b)  a  unidirectional  C-C  con^osite 
made  with  T300  fibers  and  pitch-derived  graphite  matrix;  and 
(c)  a  two-dimensional  composite  made  from  satin  weave  PlOO  fibers 
in  a  pitch-derived  graphitic  matrix. 
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As  previously  mentioned,  the  dilatometer  can  also  be  run  in  the  three- 
point-bend  mode.  In  this  arrangement,  the  sample  is  in  the  form  of  a  bar 
0.8  X  2.k  <  5  mm  (length).  The  pushrod  is  placed  on  the  center  of  the  span 
to  record  the  inelastic  deflection  as  a  function  of  temperature.  Those 
elastic  deflections  that  occur  are  of  the  order  of  a  few  microns  and  are 
neglected.  This  type  of  test  is  somewhat  uncommon;  therefore,  it  is 
appropriate  to  indicate  the  general  type  of  data  yield. 

The  inelastic  strain  of  samples  of  pyrolyzing  PAA,  phenolic,  and 
polyfurfural  alcohol  during  heat  treatment  is  shown  in  Fig.  11.  The  strain 
behavior  reveals  that  there  are  two  temperature  ranges  in  which  significant 
amounts  of  strain  energy  are  absorbed  by  the  samples.  The  first  occurs  at 
low  temperatures,  essentially  from  the  cure  temperature  to  about  1000®C, 
and  relates  to  the  formation  of  the  polymeric  carbon  submicrostructure. 

The  second  range  occurs  above  about  2000“C,  where  graphitization  is  ini¬ 
tially  observed  in  the  (formerly  glassy)  matrices  of  C-C  composites,  up  to 
the  temperature  limit  of  the  dilatometer.  Macroscopic  work  absorption  in 
this  range  is  probably  associated  with  creep  and  grain-boundary  sliding 
processes.  Such  observations  are  important  in  research  on  carboniza¬ 

tion  and  graphitization  processes.  For  further  information,  the  reader  is 
referred  to  the  literature. 
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Fig.  11.  Data  from  three-point-bend  samples  of  (a)  PAA,  (b)  phenolic,  and 
(c)  furfural  resin  samples  that  were  precured  to  350*C.  The  max¬ 
imum  shear  and  outer-fiber  tensile  stresses  were  3.5  KPa  and  0.25 
HPa,  respectively.  The  data  indicate  that  the  three  carbons  dif¬ 
fer  in  the  amount  and  temperature  range  of  strain-energy  absorp¬ 
tion.  Phenolic  resin  absorbs  more  strain  energy  during  the  early 
stages  of  pyrolysis;  PAA  absorbs  more  at  graph itizat ion  tempera¬ 
tures  (>2000*0;  and  furfural  resin  carbon  absorbs  less  at  all 
temperatures . 
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V.  CONCLUSIONS 


The  dllatometer  system  described  here  can  measure  the  thermally 
induced  strain  (deflection)  of  samples  to  a  precision  of  better  than  s5  urn 
from  room  temperature  to  3000®C.  This  precision  is  achieved  by  character¬ 
izing  the  instrumental  response  of  the  system  for  a  known  standard  that  can 
be  referenced  to  NBS  (NIST)  data,  then  algebraically  subtracting  that 
functionality  from  the  strain  of  the  unknown  sample.  The  dllatometer  thus 
emulates  a  virtual  two-pushrod  dllatometer  system.  Using  computer  data 
acquisition  and  analytical  methodology  makes  quantifying  a  sample's 
thermodynamic  responses  a  very  rapid  process. 

Non thermodynamic  properties,  such  as  the  strains  occurring  in  carbon¬ 
aceous  materials  during  their  pyrolysis,  can  also  be  measured.  Differences 
in  the  pyrolysis  behavior  of  carbon  precursors  have  been  observed.  It  is 
proposed  that  characterizing  and  analyzing  such  differences  may  shed  light 
on  the  behavior  of  precursors  in  the  matrix  of  C-C  composites  and  the 
effects  of  the  precursors  on  the  bulk  thermal  strains  of  such  composites. 
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APPENDIX  A. 


DILATOMETER  DATA  ACQUISITION  PROGRAM 

The  following  BASIC  language  code  runs  on  a  multitasking  Hewlett- 
Packard  personal  computer  under  a  BASIC  Shell  in  the  UNIX  enviroranent.  It 
is  used  for  laboratory  data  acquisition. 
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10  !««4r*****«****«*«**«*«*«**««****4>*«**«*«***«**«i**#**4>**«*** 

20  !  •* 

30  !  ASTBO  Furnace  data  acquisition  profraa  ** 

40  !  *» 

50  !  Jun  1$  1087  (Vers  2.1)  «• 

60  !  •• 

70  !  •• 

80  !  Makes  a  daughter  file  in  /usr/ASTRO  •• 

90  !  which  can  be  read  by  unix  for  VAX  transfer  ** 

100  !  •• 

110  ! 

120  SHORT  ninTiaaxTiinterTitisMefhrttthickiSthiektSero  0  INTBGBR  nopoints 
0  pointsl 

130  CLEAR  0  ALPHA  0,0  0  DIM  sa^)lt[80] 

140  D^  0  DISP  0  DISP  6  DISP 

ISO  DISP  •  Astro  Furnace  Data  progras*  0  MASS  STORAGE  IS  ■/usr/ASTRO* 

160  DISP  0  ASSIGN  16  TO  ■astro! .dat.tap* 

170  DISP  ■  calibrate  LVDT  [y/n] ■ 

180  INPUT  yesnoS  0  PRINTER  IS  16 

100  IF  yesno8=*n*  OR  yesnoS=*N*  TEEN  GOTO  220 

200  GOSUB  000 

210  GOTO  350 

240  GOTO  330 

330  DISP  6  DISP  0  DISP  ■  enter  cal  factor  (V/ua)* 

340  INPUT  calfac  0  serosO 

350  DISP  0  DISP  0  DISP 

360  DISP  ■  Enter  heating  rate  (C/sin.)* 

370  INPUT  hrt 

380  DISP  Q  DISP  ■  Enter  Pile  naae* 

300  INPUT  naaeS 
400  sla4iEN(naael) 

410  dateeSsDATEI 

420  DISP  0  DISP  ■  Saaple  I.D.  ■ 

430  INPUT  saaplS 

440  DISP  Q  DI^  ■  Enter  saaple  thickness  (aa)* 

450  INPUT  sthick 

460  DISP  0  DISP  ■  Enter  ain,aax  temperatures  (C)* 

470  INPUT  ainT,aaxT 

480  DISP  0  DISP  ■  Enter  saaple  interral  (seconds)* 

490  INPUT  interr 

500  DISP  0  DISP  8  DISP  ■  Length  of  run  (ainutes)* 

510  INPUT  tiaee 

520  tiaeees(aaxT-ainT)/hrt  !  tiae  length  of  run  IN! 

530  DISP  0  DISP  ■  Tiae  to  aaxiaua  tem>:  *;tiaeee  0  DISP 

540  nopntssINT  (tiaee/ (interr/SO))  I  total  f  of  data  points  in  run  ! 

560  GOSUB  890  !**********««***«e*«i*«**«**»*««***««start  file*************! 
570  DISP  0  DISP  0  DISP  Q  DISP  *  to  start  run,  Press  fl* 

580  ON  KEY#  !,■  start*  GOTO  600 
500  GOTO  590 

600  OFF  KEY#  1  0  DISP  0  DISP  ■starting*  0  DISP  0  DISP 
8  DISP  *  Point  tiae  teap  LVDT* 

610  ON  KEY#  2,*  HOLD*  GOTO  1380 

620  tiaasTIMB  !  starting  seconds  since  aidnight  ! 

630  dattsDATE>87000  !  starting  #  of  days  since  start  of  year 
640  !!!!!!start  of  data  loop!! 

650  ! 

660  ON  ERROR  GOTO  690 
670  ASSIGN  7  TO  *hpib.b2* 

680  OFF  ERROR  8  GOTO  695 
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12flO  OUTPUT  709  ;  'ACS*  ! 

1300  ENTBEl  709  ;  cl,c2.c3,e4,c5,c6.e7.c8,c9  ! 

1310  calf»(LL}s((el‘fe2-«’C3'»c4'»c5-fc6'»e7'»e8'fc9)/9-sero} /thick  !  T/oa 
131S  WAIT  1000  ! 

1320  NEXT  LL  ! 

1330  cslfac»(ealfm(l}>ealfa(2)-^calfa(3))/3  ! 

1340  ASSIGN  7  TO  ! 

1350  DISP  0  DISP  ■  Cal.  faetor«*;calfac  9  lEl’SN  ! 

1380  ! - 

1370  ! -  - 

1380  SISP  9  DISP  9  DISP  ■  ...Lwlding,  «nd«f3;  a>stittue<»f4* 
1390  ON  KEY#  3,*  end*  GOTO  860  9  ON  KEY#  4,*  cost.*  GOTO  1410 
1400  GOTO  1400 

1410  (VP  KEY#  3  9  OFF  KEY#  4  9  GOTO  840 
1420  recoTcr: 

1430  DISP  *  A  Tiaeottt  ha«  Occurred  on  701 . . .  RecoTering* 

1440  ASSIGN  7  TO 
1460  GOTO  650 


#### 

irvirw 

««## 

«### 

VYVf 

#### 

«### 

###« 

>###« 

-§m 
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APPENDIX  B. 


DILATOMETER  DATA  ANALYSIS  PROGRAM 

n 

The  following  IDL  language  code  is  run  on  a  VAX  11/785  computer  to 
analyze  the  data  collected  in  the  laboratory.  Data  are  transferred  to  the 
VAX  via  Kermit  file-transfer  software  (public  domain). 


n 

Interactive  Data  Language,  RSI,  Inc.,  2021  Albion  St.,  Denver,  CO. 
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SET_PLOT,0 

VTIOO 

!FANCy=0 


Version  S.OOpyrostd  baseline  7th  polynoaial 
JAN  10,  1989  with  SMOOTEING 
and  SCATTER  PLOT  and  AUTOREAD 
and  DERIVATIVB  and  LBGANDS 


PVAR=” 

print, ’  ASTRO  FILE  DATA  PLOTTER:  OVERFLOTIED  FILES  WILL  BE  OVERLAYED  IN  THE’ 
PRINT, ’  ORDER  GIVEN.  SCALE  WILL  BE  THAT  OF  TEE  FIRST  FILE  GIVEN. > 

PRINT,’  Be  sure  that  the  terainal  is  in  CAPITAL  letter  aode! 

PRINT,’  ’ 

PRINT,’  with  SMOOTHING  and  SCATTER  PI^’ 

PRINT,’  ’ 

NUMP=10 

PRINT , ’  Enter  #  oi  points  to  be  averaged  for  saoothing ’ 

read,NUMP 

PRINT,’  ’ 

PRINT, ’  NUMBER  OF  FILES  TO  BE  PLOTTED  ?  «8) ’ 

;  NUMBER  OF  FILES 

Zsl  :  A  FUG 

READ,AB 

C0DE=:1 


0RD«2 
PRINT,’  ’ 

PRINT,’  Plot  a  poljrnonial  fit  to  the  data?  [Y/N] ’ 

read,PVAR 

PRINT,’  ’ 

IP  PVAR  E«  ’Y’  THEN  BEGIN 
PRINT,’  HEAT  POLYNOMIAL  ORDER  ?’ 


READ,0RD 

ENDIF 

PRINT,’  ’ 

PRINT, ’  X-AXIS 

PRINT,’ 

Y-AXIS 

CODE  ’ 

t 

PRINT,’ 

Temperature 

Expansion 

1’ 

print, ’ 

Time 

Expansion 

2’ 

PRINT,’ 

Time 

Temperature 

3’ 

PRINT,’ 

Temperature 

Derivative  (Temp.) 

4 

(for  dilation) ’ 

PRINT,’ 

Time 

Derivative  (Time) 

5 

(for  dilation) ’ 

PRINT,’ 

Temperature 

Deflection 

6 

(for  flexure) ’ 

PRINT,’ 

Temperature 

Deflection 

7 

(for  compression) 

PRINT,’ 

Temperature 

Derivative 

8 

(for  compression) 

PRINT,’ 

Temperature 

Derivative 

0 

(for  flexure) * 

PRINT,’  ’ 

PRINT,’  Enter  CODE’ 
print,’  ’ 

READ, CODE 


NUM»INDGEN(AB) 

;  VECTOR  OF  POINT  NUMBERS 

NAME=STRARR(20,AB) 

NAMEB=” 

;  CORRESPONDING  VECTOR  OF  NAMES 

YSNO*” 

READ, ’Display  the  Directory  ’,TSN0 
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IF  YSNO  SQ  ’Y’  THEN  SPAIN, ’DIR* 


NUUEB:^ 

FOR  Isl,Afi  DO  BEGIN  ;  get  naaes  of  file*  to  be  plotted 

PRINT,  ’ENTER  NAME  OF  FILE  ’ 

READ,NAMEE 

NAME(I-1}«NA11EE 

endfor 

STDl* [-4 . 3S470B-15, 4 . 7716flE-12, -1 . 9023SE-9 ,3. 18481B-7 , -1 . 4405SE-5 , -1 . 18S9<»-3 , - . 44418 ,1.34 

;  BASELINE  FIT  DATA  JUN.  25.  1688 
STD2» [-3 . 95873B-18 , 3 . 98258B-12 , -1 . S6485B-8, 3 . 14101B-S , -1 . 36172B-2 , . 31399] 

;  COEFF’S  FOR  GRAPHITB  IN  BASAL  PLANE 
STD3» [6 . 86491E-15 , -5. 51405E-11 , 1 . 89794E-7 , -2 . 40720B-4 , 0 . 13016] 

;  BASELINE  FOR  COMPRESSION  SAMPLE 

for  i=l,ab  do  begin  ;  do  plotting... 

0PENR,1,NAME(I-1) 

a=strarr;30,i} 

B:=STRARR(60,1) 

C=fltarr{4) 

D=” 

&:DOUBLB(FLTARR(4,2000))  ;  MAKE  DATA  VECTOR  DOUBLE  PRECISION 

ON_IOERRQR, POINTS 
READF  1  A  B  G  D  E 

POINTS: ’print! (MIN (fHBRB(B  BQ  0)))/4,’  POINTS  READ’ 

ON  IOERROR,NULL 

PRINT,’  Nuaber  of  Data  Points  to  be  Plotted  ?’ 

R£AO,NUMEE 

NUM(I-1)-NUMEE 

aosE,i 


X=B(2,0!NUM(I-1)) 

Y=E(1,0:NUM(I-1)) 

XT=B(3,0:NUM(I-1)) 


DEFINE  NEW  VARIABLES  AND  SET  UP  PLOT 

MICRONS 

TIME 


PYR0FLX=(-8.6B-3)*X 


;  FOE  FLEX  SAMPLES 


Y*Y-(Y(0)+Y(l)+Y(2))/3 


SET  Y(B.T.)»0 


X*X/8 

Y1»STD1 (0) •X‘7+STD1 (1) •X‘8+STD1 (2) ♦X*5+STD1 (3) •X*4+STD1 (4) •X'3+STDl (5) *X*2+STD1 (8) •X+STD 

;  MAKE  BASELINE  DATA 

X>X*8 

Y2=STD2  (0)  •X‘6+STD2  (1 )  *X*4*STD2  (2)  •X*3+STD2  (3)  *X‘2+STD2  (4)  •I+STD2  (5) 

;  MAKE  BASAL  PLANE  DATA 
;  FOR  GR  STANDARD 


IF  (CODE  El)  7}  OR  (CODE  EQ  8}  WSH  BEGIN 

Y4=STD3  (0)  *X*5fSTD3  (1)  al'AfSIDS  (2)  «X‘34-STD3  (3)  •I*2^STD3  (4)  *X 

i  MAKE  BASELINE  FOR  COMPRESSION 

Ba>iF 


]GRID»2 

IXTITLBc’Teaperatttre  (C)  ’ 
!YTITLE*’D«lta  L/Lo  (aicrons/ca)  ’ 


U2 


!MTITL£= ’THERMAL  EZPANSION  He  fttMsphere’ 
IF  CODE  EQ  2  THEN  imTL&r’TIMB  (ainutes)  ’ 
IF  CODE  EQ  3  THEN  BEGIN 
imTLE^’TDlB  (ainutes)  ’ 

! xTITLEs ’TEMPERATURE  (C) ’ 

ENDIF 

IF  CODE  EQ  4  THEN  BEGIN 

! YTITLE= ’ dS/dT  (ua/aa/deg . C) ’ 
iMTITLEs’HIGH  TB(PERATURE  FLEIUBE’ 

ENDIF 

IF  CODE  EQ  0  THEN  BEGIN 

!MnTLE=’HIGH  TEMPERATURE  FLEXURE’ 
SYTITLE=’dS/dT  (ua/aa/C) ’ 

ENDIF 

IF  CODE  EQ  8  THEN  BEGIN 

!MTITLE»’BIGB  TEMPERATURE  COMPRESSION’ 
!yTITLE= ’ dL/dT  (ua/aa/C) ’ 

ENDIF 

IF  CODE  EQ  5  THEN  BEGIN 
!yTITLE=’dS/dt  (ua/ca/ain.)’ 
!XTITL&=’TiBe  (ainutes) ’ 

!MriTLB=’HI6H  TEMPERATURE  FLEXURE’ 

ENDIF 

IF  CODE  EQ  6  THEN  BEGIN 
!YTITLE= ’Deflection  (ua) ’ 
iXTITLEs'Teaperature  (C) ’ 

!MTITLE=’BIGH  TEMPERATURE  FLEXURE’ 

ENDIF 

IF  CODE  EQ  7  THEN  BEGIN 

!MTITLE«’HIGB  TEMPERATURE  COMPRESSION’ 
!YTITLE=’CoBpi«88iTe  Strain  (10*-3)’ 
ENDIF 


IF  (CODE  EQ  1)  or  (code  eq  4)  TEB<  BEGIN 
y»(Y-Yl+Y2)/(C(l)*10)  ;  CORRECT  DATA  FOR  DILATION  SAMPLES 

ENDIF 


IF  (CODE  EQ  7)  OS 
Y=a-Y4)/C(l) 

ENDIF 


(CODE  EQ  8)  TEEN  BEGIN 

;  SUBTRACT  COMPRESSION  BASELINE 
;  AND  MAKE  INTO  ua/aa 


IF  (CODE  EQ  6)  OR  (CODE  EQ  9)  TEEN  BEGIN 
Y»Y-PYROFLX  ;  SUBTRACT  FLEXURE  BASELINE 

ENDIF 


SBT_VIEWPORT, .11, .80, .22, .9 


IF  Z  GT  1  THEN  BEGIN  ;  USE  ’PLOT’ 

!LINETYPB=I-1 

IF  (CODE  EQ  1)  or  (CODE  EQ  6)  OR  (CODE  EQ  7) 

IF  CODE  EQ  2  THEN  OPLOT, XT, SMOOTH (Y,NUMP) 

IF  CODE  EQ  3  THEN  OPLOT,XT,X 
IF  (CODE  EQ  4)  OR  (CODE  EQ  8)  OR  (CODE  EQ  9) 
Dy=DERIV(X, (SMOOTH (Y,NUMP))) 

Dy(0)»Dy(3) 

Dy(l).Dy(0) 

OPLOT,X,Dy 


ONLY  IF  1ST  TIME  AROUND 

THEN  OFLOT,I,SMOOTH(Y,NUMP) 
;  ELSE  USB  ’OPLOT’ 

:  TIME;EXPANSION 

THEN  BEGIN 
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ENDIF 

IF  CODS  EQ  5  TEEN  BEGIN 
l>r*DEBIV(Xr,  (SllOOTB(Y,NUMP))) 

D3r(0)=Djr(3) 

l>y(l).l>y(0) 

opLOT,rr,Dr 

ENDIF 

ENDIF  ELSE  BEGIN 

!PSYM«3  ;  DO  DOTS 

VTIOO 

!LINETYPE=0 

PL0T,X,Y 

IF  (CODE  Eq  1)  OR  (CODE  EQ  6)  OR  (CODE  EQ  7)  THEN  PLOT,X,Y 
IF  CODE  EQ  2  THEN  PLOT, XT, Y 
IF  CODE  Sq  3  THEN  FLOT,XI,X 

IF  (CODE  EQ  4}  OR  (CODE  EQ  8)  OR  (CODE  BQ  8)  THEN  BEGIN 
D7sDERIV(X,  (SMOOTH  (Y,NU|(P) )  ] 

Dy(0)=Dy(3)  ;  GENERAL  DERIVATIVB  CURVES 

Dy(l)-Dy(0) 

PLOT,X,l>y 

ENDIF 

IF  CODE  EQ  5  TEEN  BEGIN 

Dy=DERIV(rT, (SMOOTH (Y,NUMP))) 

Dy(0)=Dy(3) 

Dy(l)=Dy(0) 

PLOT, XT, Dy 

ENDIF 

!PSYM=0 

IF  PVAR  NE  ’Y*  TEEN  BEGIN  ;OPLOT  SMOOTHED  CURVE 

SET  XY 
!FSTM=0 
!LIHBTYPB~0 

ip  (CODE  EQ  1)  OR  (CODE  EQ  8)  OR  (CODE  BQ  7}  TIEN  OPLOT, X. SMOOTH a, NUMP) 
IF  CODE  EQ  2  THEN  OPLOT, XT, SMOOTH (Y, NUMP) 

IF  CODE  EQ  3  THEN  OPLOT, XT,X 

IF  (CODE  EQ  4)  OR  (CODE  EQ  8)  OR  (CODE  BQ  6)  TEEN  BEGIN 
Dr=DBEIV(X,  (SMOOTH  nr,  NUMP) ) ) 

Dy(0)»Dy(3) 

Dy(l)=D7(0) 

OPLOT, X,Dy 
ENDIF 

IF  CODE  EQ  5  THEN  BEGIN 

Dy=DERIV(XT, (SMOOTH (Y, NUMP))) 

Dy(0)»Dy(3) 

Dy(l)«Dy(0) 

OPLOT, n,Dy 
ENDIF 

ENDIF 

IF  PVAR  EQ  ’Y’  THEM  BEGIN  ;  OPLOT  POLYNOCUL  CURVE  IF 

COBFF-POLY  FIT(X,Y,ORD,GOG)  ;  PVAR»Y 
IPSYM-0  - 
OPLQT,X,6GG 
ENDIF 


ENDELSE 

BNDFOR 


!NOEEAS«l 

SET  ¥IBiP0ST,O,l,O,l 
SET  Ilf, 0,1,0, 1 
ABCl^.82 

XY0UTS,.85,.88,’Naaes:> 

FOR  LUl,Afi  DO  BEGIN 
NNNsFINDGEN(AB>l} 

!LIHETTPE=NNN(LL-1) 

OPLOT, [ . 82 , . 86] , [ABCD- (LL-1) * .05 , ABO)- (LL-l) * .05] 
XYOUTS, .875, ABCD- {LL-l)*.05-.01,NAiE(LL-l) 

ENDFOR 
SET  XY 
SBT'VIBiPORT 
!NOERAS=0 

SET  VIEWP0ET,O, 1,0,1 
SBT“XY,0,1,0,1 

XY0UTS,.17,.l,’N««e,  Run  Date;  ' 
rrOUTS,.39,.l,A 

XYOUTS, .17, .05, ’Heating  Rate  (C/ain.):  ’ 

XYOUTS, .455, .05, C(0) 

XYOUTS,. 7,. 05, ’Lo  (ob) ;  ’ 

XYOUTS,. 79,. 05, C(l) 

XYOUTS, .17, .075, ’Saaple  ID:  ’ 

XYOUTS, .37, .075, B 

SET  XY 

SET”VIB1P0BT ,.11, .80, .22, .9 


5  *♦♦**♦•*♦*•♦♦*••♦*•••**••*•*•**•*♦*♦♦•*♦*♦*•*•••***•*•**••*••••**•• 

iPLOTTING  SECnON 


PRINT, ’  SAVE  THIS  ?  [Y/N] ’ 

YN=” 

EEAD,YN 

IF  YN  Eq  ’Y’  THEN  BEGIN 
Z=1 

for  i=l,ab  do  begin 
0PENR,1,NAMB(I-1) 

At=STEABR(30,l} 

BsSTRARB(60,1} 

C*fltarr(4) 

D=” 

E^DOUBLE  (FLTARR  (4 ,  NUli  (I-l)  -^1}  ) 

READF,1,A,B,C,D,E 

aOSE.l 

I=E(2,0:NUM(I-1)) 


:  IF  HARD  COPY  IS  DESIRED 
;  do  plotting... 


;  MAKE  DATA  VECTOR  DOUBLE  PRECISION 

;  DEFINE  NEW  VARIABLES  AND  SET  UP  PLOT 


Y*B(1,0:N0M(I-1)) 

XT=B(3,0:N01I(I-1)) 

mOPLX-(-6.8B-3)*X 


;  MICRONS 
;  TIMB 


;  POE  FLEX  SAMPLES 
Y«Y-(Y(0)+Y(l)+Y(2))/3  ;  SET  ZERO  OP  CURVE 


X=X/8 

Yl-STDl (0) •X*7*STD1 (1) ♦X*6+STD1 (2) •X*5>STOl (3) *X“4+STD1 (4) •X‘3+STD1 (5) •X‘2+STD1  ^5) •X^STDl ( 

;  MAKE  BASELINB  DATA 

X=X*8 

Y2=STD2 (0) ♦X*5+STD2 (1) *X‘4+STD2 (2) •X*3+STD2 (3) •X‘2+STD2 (4) *X>STD2 (S) 

;  MAKE  BASAL  PLANE  DATA 
:  POR  GR  STANDARD 


IF  (CODE  EQ  1}  or  (code  eq  4)  THEN  BEGIN 
Y»(Y-Yl+Y2)/(C(i)*10)  ;  CORRECT  DATA  FOR  DILAHON  SAMPLES 

ENDIF 


IF  (CODE  E«  7)  OR  (CODE  EQ  8)  THEN  BEGIN 
Y«(Y'Y4)/C(1)  ;  SUBTRACT  COMPRESSION  BASELINB 

;  AMD  MAKE  INTO  rm/wm 


ENDIF 


IF  (CODE  EQ  6)  OR  (CODE  EQ  9)  TEEN  BEGIN 
Y=Y-FYROFU  ;  SUBTRACT  FLEXURE  BASHINB 

ENDIF 


Ifeacy  »  2 
set_plot,4 

SET^VIBiPOBT, .11, .80, .22, .9 

IF  Z  6T  1  THEN  FEGIN  ;  USE  ’PLOT’  ONLY  IF  1ST  TIME  AROUND 

!LINBTYPE=I-1 

IF  (CODE  EQ  1)  or  (CODE  EQ  6)  OR  (CODE  EQ  7}  THEN  QPLOT,X, SMOOTH (Y,NUMP) 

;  ELSE  USE  ’OFLOT’ 

IF  CODE  EQ  2  THEN  OFLOT,XT,SMOOTHa.NUMP)  ;  TIME: EXPANSION 

IF  CODE  EQ  3  THEN  OPLOT,XT,X 

IF  (CODE  EQ  4)  OR  (CODE  EQ  8)  OR  (CODE  EQ  9}  THEN  BEGIN 
Djr-DERIV(X,  (SMOOTH  (Y,NUMP))) 

Djr(0)»Dy(3) 

Djf(l)=Dy(0) 

0PL0T,X,D7 

ENDIF 

IF  CODE  EQ  5  THEN  BEGIN 
Dy«DBRIV(XT, (SMOOTH (Y,NUMP))) 

D3r(0)=Dy(3) 

D3r(l)»Dr(0) 

0FL0T,XT,D7 

ENDIF 

ENDIF  EL  SE  BEGIN 

!PSYlt:3  ;  DO  DOTS 

VTIOO 

ILINETYFE-O 

PLOT.X.Y 

IF  (CODE  EQ  1)  OR  (CODE  EQ  6}  OR  (CODE  EQ  7)  THEN  FLOT,X,Y 


il6 


IF  CODE  B9  2  THEN  PLOT,XT,Y 

IF  CODE  EQ  3  THB«  PLOT, XT, X 

IF  (CODE  M  4)  OR  (CODE  EQ  8)  OR  (CODE  EQ  8)  THEN  BEGIN 
Dy>^RIV(X,  (SMOOTH (Y,Nl»F)}) 

D7(0)a4>y(3}  ;  GENERAL  DERIVATIVE  CURVES 

Dy(l).Dy(0) 

PLOT,X,Dy 

ENDIF 

IF  CODE  EQ  5  TH^  BEGIN 

Dr=DERIV(XT,  (^OOTE(Y,NUMP})} 

Dy(0)=Dy(3) 

Dr(l)*Dy(0) 

PLOT,XT,Dy 

ENDIF 
•PSYIM) 

IF  PVAR  NE  *Y*  THEN  BEGIN 
SET  XY 
»PSYM=0 
!LINBTYPE=0 

IF  (CODE  EQ  1)  OB  (CODE  EQ  6)  OR  (CODE  EQ  7) 

IF  CODE  EQ  2  THEN  OPLOT, XT, SMOOTH (Y,NUMP) 

IF  CODE  EQ  3  THEN  OPLOT,XT,X 
IF  (CODE  4}  OR  (CODE  EQ  8)  OR  (CODE  EQ  8} 

Dy«DEEIV(X, (SMOOTH(Y,NUMP))) 

Dy(0)=Dy(3) 

D7(l)=Dy(0) 

OPLOT,X,Dy 
ENDIF 

IF  CODE  EQ  5  TEEN  BEGIN 

Dy=DERIV(XT, (SMOOTH(Y,NUMP))) 

Dy(0)*Dy(3) 

Dy(l)=Dr(0) 

0PL0T,XT,Dt 
ENDIF 

ENDIF 

XYOUTS, iCXMAX*0.843, <CYMAX*1.02, !STIME 

IF  PVAR  EQ  ’Y’  THEN  BEGIN  ;  OPLOT  POLYNOMIAL  CURVE  IF 

COBFF»POLY  FIT(X,Y,ORD,GGG)  ;  PVAR*Y 
!PSYM=0  ~ 

OPLOT, X,GG6 
ENDIF 

Z=2 

SNDELSE 

SNDFOR 

endif 

IF  CODE  EQ  1  THEN  BEGIN 
IF  Z  EQ  2  THEN  BEGIN 
IF  AB  EQ  1  THEN  BEGIN 

PRINT,  ’  TO  CALCULA1X  ALPHA,  ENim  NUMBER  OF  REGIONS  TO  BE  FITTED.  ’ 


; OPLOT  SMOOTHED  CURVE 

TEEN  OPLOT,X.SHOOTB(Y,NUMP) 
THEN  BEGIN 
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PBIMT,’  TO  END,  ENTER  ZERO' 

ANSVEEUO 

READ.ANSIER 

if  ANSWER  6T  0  TEEN  BEGIN 
FOR  I«l, ANSWER  DO  BEGIN 
PRINT,'  ' 

PRINT, ’  ENTER  TEMP  RANGE  (LOWBR,UPPER) ' 

L0WER«5 

OPPERsS 

READ, LOWER, UPPER 

LOWBR»IC[N(WEERE(Z  GT  LOWER))  ;  GET  SUBSCRIPTS  OP  TEUP  RAN(3 

UPPERaMAX(WHEBB(Z  LT  UPPER)) 

ZFIT»X(LOWSt:UPPER)  ;  DEFINE  NEW  DATA  VECTORS  FOR 

YFIT-Y(LOWER:UPPER)  ;  ALPHA  FIT 

COEPF*POLY  FIT(XPIT,YPIT,1,LINE)  ;  FIT  LINE  TO  DATA 

iLINETYPSa 

OPLOT,IFIT,LINE 

ALPHA«COEFF(1)*100  ;  MAKE  ALPHA  INTO  IdCRONS  PER  METER 

XyOUTS,XCUPPER)  ,Y(UPPER) , ALPHA 

SBT_PLOT,0 

SBTjnBWPORT, .11, .80, .22, .8 

!FANCY=0 

OPLOT,XFn,LINE 

XYOUTS,X(UPPER)  ,Y(UPPER)  .ALPHA 
SET_PL0T,4 

SBTjnffiWPOBT, .11, .80, .22, .8 

rFANCY=3 

ENDFOR 

ENDIF 

ENDIF 

ENDIF 

ENDIF 

IF  YN  EQ  'Y*  THEN  BEGIN 
INOERASsl 

SET  VIBWP0RT,O, 1,0,1 

SET''Xy,0,l,0,l 

ABCE-.82 

XYOUTS,  .85,  .86, 'Nuns; ' 

FOR  LLsl.AB  DO  BEGIN 
NNN>FINDGBN(AB*1) 

!LINETYFB:4QiN(LL-l) 

OPLOT,  [.80,  .87] ,  [ABCD-(LL-1)*.06,ABCD-{LL-1)*.06] 

XYOUTS, .875, ABCD-(LL-1} • .05-. 01 , NAME (LL-1) 

ENDFOR 

XYOUTS,  .17,  .1, 'Nue,  Run  D»te:  ’ 

XYOUTS,. 31,. 1, A 

XYOUTS, .17, .05, 'Heating  Rat«  (C/ain.):  * 

XYOUTS, .33, ,05,0(0) 

XYOUTS,. 6,. 06, 'Lo  (c«);  ’ 


XyOUTS,.67,.06,C(l) 
XyOUTS,.l7,.07S,'Sa«ple  ID:  ’ 
IY0UTS,.29,.075,B 

'SBTJW 

SET  VIEirP0RT,.ll,.8O,.22,.6 
!N0^AS»O 

VTIOO 

PRINT, *  LASER  PRINTING  TEE  FILE  NOW. . 
WAIT, 2 
SET  PLOT,0 
VTlSo 

SPAWN,  *LASE  DOWN  (}tlSPLOT.LIS' 

ENDIP 


TECHNOLOGY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for  national  security 
programs,  specializing  in  advanced  military  space  systems.  The  Corporation's  Technology 
Operations  supports  the  effective  and  timely  development  and  operation  of  national  security 
systems  through  scientific  research  and  the  application  of  advanced  technology.  Vital  to  the 
success  of  the  Corporation  is  the  technical  st^fs  wide-ranging  expertise  and  its  ability  to  stay 
abreast  of  new  technological  developments  and  program  support  issues  associated  with  rapidly 
evolving  space  systems.  Contributing  capabilities  are  provided  by  these  individual  Technology 
Centers: 

Electronics  Technology  Center:  Microelectronics,  solid-state  device  physics, 

VLSI  reliability,  compound  semiconductors,  radiation  hardening,  data  storage 
technologies,  infrared  detector  devices  and  testing;  electro-optics,  quantum  electrcmics, 
solid-state  lasers,  optical  propagation  and  communications;  cw  and  pulsed  chemical 
laser  development,  optic^  resonators,  beam  control,  atmospheric  propagation,  and 
laser  effects  and  countermeasures;  atomic  frequency  standards,  applied  laser 
spectroscopy,  laser  chemistry,  laser  optoelectronics,  phase  conjugation  and  coherent 
imaging,  solar  cell  physics,  battery  electrochemistry,  battery  testing  and  evaluation. 

Mechanics  and  Materials  Technology  Center:  Evaluation  and 
characterization  of  new  materials:  metals,  alloys,  ceramics,  polymers  and  their 
composites,  and  new  forms  of  carbon;  development  and  analysis  of  thin  Olms  and 
deposition  techniques;  nondestructive  evaluation,  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  development  and  evaluation  of 
hardened  components;  analysis  and  evaluation  of  materials  at  cryogenic  and  elevated 
temperatures;  launch  vehicle  and  reentry  fluid  mechanics,  heat  transfer  and  flight 
dynamics;  chemical  and  electric  propulsion;  spacecraft  structural  mechanics, 
spacecraft  survivability  and  vulnerability  assessment;  contamination,  thermal  and 
structural  control;  high  temperature  thermomechanics,  gas  kinetics  and  radiaticui; 
lubrication  and  surface  phenomena. 

Space  and  Environment  Technology  Center:  Magnetospheric,  auroral  and 
cosmic  ray  physics,  wave-particle  interactions,  magnetospheric  plasma  waves; 
atmospheric  and  ionospheric  physics,  density  and  composition  of  the  upper 
atmosphere,  remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared 
astronomy,  infirared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere;  effects 
of  electranagnedc  and  particulate  radiations  on  space  systems;  space  instnunentadon; 
propellant  chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
atmospheric  chemical  reactions,  atmospheric  optics,  light  scattering,  state-speciflc 
chemical  reactions  and  radiative  signatures  of  missile  plumes,  and  sensor  out-of-field- 
of-view  rejection. 


